Exact solutions of an elliptic Calogero— Sutherland model 
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A model describing N particles on a line interacting pairwise via an elliptic function potential in 
the presence of an external field is partially solved in the quantum case in a totally algebraic way. 
As an example, the ground state and the lowest excitations are calculated explicitly for N = 2. 

PACS numbers: 03.65.Fd, 71.10.Pm, ll.10.Lm 



It is well known that the class of exactly solvable prob- 
lems does not include most physical problems. The devel- 
opment of computer science in the last decades has made 
possible the use of numerical methods to approximate 
exact solutions in a wide variety of situations. Yet, the 
study of exactly solvable models still deserves attention, 
not only because the knowledge of exact solutions can be 
used to test approximate methods, but also in its own 
right, due to the simplicity and mathematical beauty of 
the models, and the wide range of connections with other 
fields of physical and mathematical research. 

This is illustrated by the renewed interest in the 
Calogero-Sutherland (CS) models of interacting parti- 
cles in one dimension, which have been recently applied 
to many different fields such us quantum spin chains with 
long range interaction jj| , random matrix theory Q , frac- 
tional statistics and anyons ||, Yang-Mills theories j|, 
quantum Hall liquids [||, soliton theory j^], vicinal sur- 
faces in crystals jjj , and black holes || . 

The first example of a non-trivial integrable quantum 
many-body problem was found by Calogero || , and con- 
sists of a system of identical nonrelativistic particles in- 
teracting pairwise through an inverse-square potential 
v(r) = r~ 2 , so that 

N N 

H N = -J2€ k +9^v(x j -x k ). (1) 

fc=l j,k=l 

By integrable we mean here that a complete commuting 
set of constants of motion can be explicitly constructed. 
Soon afterwards, Sutherland jio) established the integra- 
bility of the model (Q) with an inverse sine-square inter- 
action v(r) = sin~ 2 (r). 

The most general interaction potential v for which 
the Hamiltonian (|]) is known to be integrable is the 
Weierstrass V function, which includes the rational and 
trigonometric cases as special limits. The integrabil- 
ity of this potential was proved in the classical case by 
Calogero and Perelomov by means of a Lax pair repre- 
sentation p"l| ] , and its explicit integration was performed 
by Krichever, [Q. Olshanetsky and Perelomov [|l3| later 
showed that all these models have an underlying alge- 
braic structure based on root systems of A n algebras, 



and that integrable models associated to other root sys- 
tems also exist. In the models treated in Ref. |D| the 
integrals of motion are related to the radial parts of the 
Laplace-Beltrami operator on a symmetric space asso- 
ciated to the given root system. These integrable mod- 
els are obtained from the projection of free motion on a 
higher-dimensional manifold. 

However, integrable Hamiltonians are not necessarily 
solvable, i.e., we might not be able to find explicitly their 
spectrum and eigenfunctions. The models with inverse- 
square and inverse sine-square interaction are known to 
be solvable, and much literature has been devoted to the 
study of their eigenfunctions |L4j, but the more general 
model with the Weierstrass V potential is considerably 
more difficult. In fact, very few explicit solutions are 
known in this case fllH , and only for a low number of 
particles. 

The purpose of this Letter is to present a model of N 
particles on a line with elliptic pairwise interaction in an 
external field for which a finite number of eigenvalues and 
eigenfunctions can be computed algebraically. We shall 
only sketch here the main ideas behind the proof of this 
result, referring the reader to our previous work |]l6| for 
a more complete description of the method used. 

Consider the iV-body quantum Hamiltonian 

N 

H N = -Y,d 2 Xk +V N (x) (2) 
fc=i 

with potential 

N N 

Vat(x) = Cm J2 ?{ x k + i0) + ib{b V{2x k ) 

k=l k=l 
N 

+ 0(0 - 1) \- V ( X 3 + X k) + V(. X 3 ~ X k)] , (3) 

],k=l 
3^k 

where a and b are positive real parameters, m is a non- 
negative integer, 

c m = 2[2b + m + a(N - 1)] 

x [26 + 2m + 2a(N - 1) + l] > , (4) 
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and V(z) = V{z\g2,g^) denotes the Weierstrass V func- 
tion with invariants 52 , 53 G R- If 32 and 33 satisfy the 
inequality g\ > 27g|, then V(z) has two fundamental 
periods 2a and 2i/3 which are real and purely imaginary, 
respectively In this case V(x + i/3) is real and regu- 
lar (analytic) for all real values of x, with real period 2a. 
On the other hand, V{x) is real for real x and diverges as 
(x — 2na)~ 2 when x tends to an integer multiple 2na of 
the real period 2a. Thus, the configuration space for the 
Hamiltonian can be taken as the bounded region 

of R N 



< xn < xjv-i < • ■ • < x\ < a 



(5) 



Since the potential (|3[) is confining in this region, the 
spectrum of iJ/v is purely discrete, and the boundary 
condition satisfied by its eigenfunctions , 0fe( x ) is their 
vanishing on the boundary of (||). 

The potential (j|) with c m = is of Cjv type Q 



N N 

46(6-1) ^ v(2x k )+a(a-l) % 



k=l 



j,k=l 



[v(Xj-X k )+v(Xj+X k )] , 



with interaction potential v(r) = V{r). The term pro- 
portional to c m in ([?]) can be viewed as the contribution 
of an external field with potential V(r + i/3). 

We shall now show that, when the parameter to is 
a non-negative integer, one can algebraically compute 
a finite number (depending on to, see @ below) of 
eigenvalues and eigenfunctions of the Hamiltonian (^)- 
(|3|). These algebraic eigenfunctions have the form 



t/>k(x) = /x(x) Xfc(z)- 



(6) 



where 



M (x) - [] [7>fo + i/3) - P(x fe + tf)]" II ^'(^ + ^ > 

(7) 

and Xfc (z) is a suitable completely symmetric polynomial 
of degree at most Nm in the variables 



(8) 



The exact solutions of the trigonometric and rational C n 
(in fact, BCn) models also assume the form (^|), but in 
this case \x can be factorized over the system of positive 
roots, i.e., it has the form 



j<k 



x k )f( Xj +x k )] a l[[f(2x k )] 1 



(9) 



where f{x) — x in the rational case and f(x) — sinx in 
the trigonometric case. Moreover, in both cases /i coin- 
cides with the ground-state wave function of the system. 
By contrast, the function [i in cannot be factorized 



over the system of positive roots for any function fix), 
and is not the ground-state wave function of the Hamil- 
tonian (||)-([3|). As a matter of fact, it was shown in Q 
that the most general potential allowing for the factor- 
ization (H)-@ does not include the elliptic case. This is 
one of the reasons why it has been so difficult to obtain 
explicit solutions of the elliptic CS models. 

When the parameters a and b are positive, the func- 
tions (||) are regular in the region (||), and they automat- 
ically vanish on its boundary on account of the identities 
V'(i/3) = V'(a + i/3) = (see Ref. |L7|). Thus, to show 
that ip k in Eq. (^) is an eigenfunction of Hn we only 
have to check that it satisfies the Schrodinger equation 
(Hn — Ek)ip k — in the open region (||). Equivalently, 
Xk must be a solution of the equation (Hn — EkjXk = 0, 
where the gauge Hamiltonian Hn is defined by 



Hn — t l Hn M ■ 



(10) 



Note that, by the standard properties of the Weierstrass 
function []l7j, fi does not vanish in the region (||). 

It can be shown that, provided that to is a non-negative 
integer, Hn preserves the finite-dimensional polynomial 
space 



M m = span jr^ 2 • ■ • : ^ h < ml 



(11) 



where 



T k 



z tl z l2 ---z lkl l<k<N, (12) 



il<i2<-<ik 



are the elementary symmetric functions of the variables 
Zfe. This is due to the fact that, when H n is written 
in terms of the symmetric variables t\, . . . ,tn, it can be 
expressed as a quadratic combination of the generators 
of sl(N + 1) in the representation 

V k = d Tk , Mjk = Tj d Tk , U k = T k ^to - ^2 T * ^ 5 



j,k = l,2,...,N. 



(13) 



Since these generators obviously preserve the subspace 
(|TT|), so does the gauge Hamiltonian Hn- It follows that 
Hn has (at most) 



dim .M m = 



N 



(14) 



eigenfunctions \k lying in M m , which can be al- 
gebraically computed, along with their correspond- 
ing eigenvalues, simply by diagonalizing the finite- 
dimensional matrix of H^\M m - The elements of M. m , 
being polynomials in the symmetric variables r k of de- 
gree at most to, are symmetric polynomials in z of degree 
not greater than Nm. Thus the original Hamiltonian 
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(H) - © possesses (at most) dim.M m algebraically com- 
putable eigenfunctions of the form (||)-(0), with Xfc( z ) a 
symmetric polynomial of degree at most Nm, as claimed. 
Note, however, that there are other eigenfunctions of (|^) 
which do not belong to this algebraic sector. An interest- 
ing open problem is to analyze the position in the spec- 
trum of the algebraic eigenvalues and their degeneracy. 
The algebra sl(N +1), which plays a fundamental role 
in the partial integrability of the Hamiltonia n (p| )-(p[), is 
sometimes called a hidden symmetry algebra foW, since in 
this case the Hamiltonian need not be a Casimir element. 

For simplicity (indeed, no conceptual difficulties arise 
for higher values of N and to), let us consider the problem 



for N — 2 and to = 2, for which the potential reads 

2 2 

V 2 (x u x 2 ) = c 2 J2 Vix k + i0) + 46(6 - 1) 7>(2^fc) 
fc=i fc=i 
+ 2a(a - l)[V(xi + x 2 ) + V{x x - x 2 )\ . (15) 

Note that this is intrinsically a two-body problem, since 
the potential ( fl5| ) is not translation invariant. The num- 
ber of algebraic eigenstates is at most dimA^2 = 6, and 
the matrix of the restriction H 2 \_m 2 with respect to the 
canonical basis {1, tj., t 2 , t\ , tit 2 , t| } of A4 2 is given by 



/ g 2 {2a + 2b + l) -2ag 3 

16a + 246 + 20 g 2 (b+ 1/2) 

8a + 246 + 12 

8a + 126 + 14 

8a + 126 + 14 

V 



4ff3 

4ff 2 (a + 6+l) 



16(a + 36 + 3) 






2gs(l-o) 
52 (2a + 26 + 5) 
52(6+1/2) 


8a + 246 + 28 






-453(0+1) 

2g 3 
32(26 + 3) 




(16) 



/ 



The eigenvalues of this matrix are the algebraic ener- 
gies of the physical Hamiltonian H 2 , and its eigenvectors 
give the components of the corresponding functions Xk 
in Eq. ([|) with respect to the canonical basis of M 2 . 

Take, for instance, the following values of the coupling 
constants and the invariants of V(x): 

a = 2, 6-3/2, 52 = 3, g 3 = 2/3 , (17) 

for which the half-periods are 

a = 1.31523, i/3 = 1.61809 i. 

The matrix (Hq) has six real eigenvalues, so that the 
Hamiltonian p^-(0) possesses six algebraic energies, 
which have been listed in Table |[ In this table the (ex- 
act) algebraic eigenvalues of H 2 have been compared with 
a numerical approximation of the six lowest energy levels 
obtained by solving the Schrodinger equation in the trian- 
gle < x 2 < x\ < a using a finite clement method. The 
approximate agreement between the two columns of Ta- 
ble I shows that the algebraic eigenvalues are in this case 





Exact 


Numerical 




-86.5484 


-86.40 


£1 


-43.2786 


-42.96 


E2 


-10.0288 


-9.64 


E3 


12.5045 


13.11 


E 4 


46.2657 


46.97 


E 5 


81.0857 


81.34 



TABLE I: Algebraic eigenvalues Ei of the Hamiltonian 
([[5]) and their numerical approximations. 



the lowest energy states of the system, although there is 
no guarantee that this should still be true in the general 
case. The polynomials Xk corresponding to the six alge- 
braic eigenfunctions of the potential (^5|)~(|l^) are given 
in Table ||. In Figs. [l]-|| we present a plot of the ground 
state and the first two excited wave functions of the sys- 
tem. 





Tl 


T2 


r? 


T\T2 


rf 


Xo 


-3.0585 


7.1643 


2.2349 


-9.778 


9.0382 


Xi 


-3.3008 


-10.917 


2.5612 


13.462 


-24.885 


X2 


0.40412 


5.1067 


-2.2912 


-0.77221 


6.1599 


X3 


-3.9273 


-37.989 


2.0457 


14.734 


94.264 


X4 


1.6422 


-2.7456 


-0.17611 


-10.472 


-18.108 


X5 


3.8942 


8.3079 


3.6675 


15.1 


14.898 



TABLE II: Polynomials \k corresponding to the algebraic 
eigenfunctions ipu = MXfc 

of the Hamiltonian (§)-©. In all 
cases, the coefficient of 1 has been normalized to unity. 

The fact that ipo = /j,Xo is the ground state of the 
system is immediately apparent if we note that xo can 
be expressed in terms of the variables z k as 

Xo = 2.2349(z 2 - 0.82835)(z 2 - 0.54017) 

- 9.778 zi [z 2 - 0.79769)(z 2 - 0.39212) 

+ 9.0382 zl(z 2 - 0.75385)(z 2 - 0.32801). (18) 

For the values of the invariants g 2 and 33 given in ([T^), 
we have (T^j 

e 3 = -0.72011 < z k = V{x k + i/3) < e 2 = -0.240851 . 
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tpo{xi,x 2 ) 




351 1 




FIG. 1: Ground-state wave function of the potential (fl5|). 



FIG. 3: Second excited state wave function of the potential 



ipi{xi,x 2 ) 




We thus see from the previous expression that \o is pos- 
itive everywhere. Since, by ([)]), fi has no zeros in the 
open triangle < xi < x\ < a, it follows that tpo does 
not vanish in this triangle. 

In conclusion, a finite number of eigenvalues and eigen- 
functions of a quantum Hamiltonian describing N parti- 
cles on a line with elliptic interaction in the presence of 
an external field have been explicitly calculated by an al- 
gebraic method independent of the usual approach based 
on root systems. 
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